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ABSTRACT 

We report rotational periods for 16 young brown dwarfs in the nearby Upper Scorpius as¬ 
sociation, based on 72 days of high-cadence, high-precision photometry from the Kepler space 
telescope’s K2 mission. The periods range from a few hours to two days (plus one outlier at 
5days), with a median just above one day, confirming that brown dwarfs, except at the very 
youngest ages, are fast rotators. Interestingly, four of the slowest rotators in our sample exhibit 
mid-infrared excess emission from disks; at least two also show signs of disk eclipses and accre¬ 
tion in the lightcurves. Comparing these new periods with those for two other young clusters 
and simple angular momentum evolution tracks, we find little or no rotational braking in brown 
dwarfs between 1-10 Myr, in contrast to low-mass stars. Our findings show that disk braking, 
while still at work, is inefficient in the substellar regime, thus provide an important constraint 
on the mass dependence of the braking mechanism. 

Subject headings: 


1. Introduction 


Rotation is a key parameter in stellar evolution, 
and is directly linked to many fundamental astro- 
physical processes, including magnetic field gener¬ 
ation, stellar winds, star-disk interaction, and bi¬ 
nary formation (Herbst et al.||2007 Bouvier et al. 


20141. Because rotation drives magnetic activity. 


it is also relevant for ’space weather’ in the stel¬ 
lar surroundings and thus for planetary habitabil¬ 
ity (Giidel 2007 Vidotto et al. 20141. The rota¬ 
tion period, derived from the periodic photometric 
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modulation induced by surface spots, is one of the 
rare stellar parameters that can be measured with 
high precision 1%) for large samples. 

Large numbers of periods have been measured 
for low-mass stars with spectral types F to M of all 


ages from ~ 1 Myr 

to several Gyrs (e.g. Herbst et 

al. 112002 

Lamm et al.||20051 |lrwin et al.||2008| |Mc- 

Quillan et al.||2014 

I, and they provide very useful 


constraints for angular momentum evolution mod¬ 


els (e.g. Gallet & Bouvier 20151. Other portions 


of the mass-age parameter space are still poorly 
explored. One example is the regime of BDs, sub- 
stellar objects with masses less than ~ 0.08 Mq. 
The only regions with significant samples of mea¬ 
sured BD periods are clusters in Orion, the ONC 


(Rodrfguez-Ledesma et al. 

2009 

1 crOri ( 

Scholz & 

Eisloffel 

2004 Gody & Hillenbrand 

2010 

1, and 

eOri ( 

Scholz & Eisl6ffel||2005 

I, at ages of 1-5Myr. 


Like stars, BDs start with a range of periods from 
hours up to 5-10 days, and spin up due to pre-main 
sequence contraction, but in contrast to stars they 
maintain fast rotation with negligible spindown for 
several Gyrs. 

Due to their intrinsic faintness, the small pho- 
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tometric amplitudes, and the wide range of peri¬ 
ods at young ages, it is challenging to obtain a 
reliable picture of the rotational evolution of BDs 
with ground-based observations alone, which fea¬ 
ture typical daytime gaps and are therefore biased 
in the period sensitivity. There are good indica¬ 
tions that disk braking is less efficient in the very 
low mass regime (e.g. Lamm et al. [2005 1. But the 
lack of substantial and unbiased period samples 
for a range of ages prohibits us from putting firm 
constraints on models. 

Here we present first period measurements for 


BDs obtained with Kepler’s K2 mission (Howell 
et al.|[MT4). In campaign 2, the K2 field included 


large parts of the Upper Scorpius (UpSco) star 
forming region (age 5-10 Myr), which harbours a 


huge population of BDs (e.g. Lodieu et al. 2008 


Slesnick et al.|[2008{ [Lodieu et al. 2011 Dawsonj 
et aLpOld ). With a cadence of 30 min and a time 
baseline of 72 d, the K2 lightcurves cover the entire 
period range for young BDs. During the course of 
its mission, K2 is expected to fill the period-age 
diagram for BDs with dozens of new datapoints; 
this study is the first step towards a new census of 
BD rotation. 


2. Data processing 


tor and the measured flux to correct for artifacts 
induced by the pointing drift of the spacecraft. 

Our data reduction proceeds as follows. First 
we used PyKE’s kepmask to define sub-image 
apertures (~ 2 — 3x2 — 3 pixels) from the 
respective target’s postage stamp image (typ¬ 
ically 10 — 12 X 10 — 12 pixels), then extract 
simple aperture photometry (SAP) with kepex- 
tract. Next, we detrend and normalize the raw 
lightcurves with kepflatten using a low-order poly¬ 
nomial (n ~ 3 — 5) and, finally, to correct for the 
motion-induced systematic effects we use kepsff. 
As an illustrative example, on Fig. we show 
a section of the SAP lightcurve (background- 
corrected) of EP1C204614722 (top panel), the two 
moment centroids (middle panels) and the corre¬ 
sponding kepsjf-coTrected lightcurve (lower panel). 
Throughout the data-reduction process we appro¬ 
priately adjust the relevant input parameters on a 
target-by-target basis. We note that we used the 
entire campaign 2 dataset, including the majority 
of cadences from the second half of the campaign 
where most of the data flags are non-zero due to 
the photometer’s local detector electronics par¬ 
ity errors, as these do not affect the quality of 
the data for the purposes of detecting rotational 
variability. 


We observed brown dwarf^ in UpSco as part 
of campaign 2 of the K2 mission. The targets 
were selected from the census by [Dawson et al.j 
(2013) based on the UKIDSS survey. Analysis of 
mid-infrared data shows that ^ 20% still harbour 
circum-sub-stellar disks. 51 of the full sample of 
116 have been covered by K2. 

We downloaded the data, in the form of long- 
cadence target pixel files, from the Mikulsi Archive 
for Space Telescopes (MAST). To extract aperture 
photometry, detrend and correct for systematic 


effects we used the latest version of PyKE (Still 
h Barcl^|2012 ). This is an open-source PyRAF 


package originally designed for custom manipula¬ 
tion of the Kepler data and now also including the 
self-flat-fielding procedure of[Vanderburg fc John- 


son 


(2014) optimized for processing data from the 
K2 mission. The procedure utilizes the relation 
between the position of the target on the detec- 


^ Some of our objects may be very low mass stars with masses 
slightly above the substellar limit; for convenience we will 
call all targets brown dwarfs. 


3. Lightcurve analysis 


All 51 lightcurves were visually inspected. 16 
of them show clear signs of periodic variability. 
In most cases the periodicity is obvious and peri¬ 
ods ranging from a few hours up to 5 d are esti¬ 
mated by eye. To measure the periods, we used 
the autocorrelation function, which has been es¬ 
tablished as a useful tool for period analysis in Ke¬ 
pler lightcurves by jMcQuillan et aL (2014). Com¬ 
pared to Fourier (Lomb-Scargle) or phase-based 
techniques (e.g., phase dispersion minimisation), 
it has the advantage that it is robust against 
phase/amplitude changes as well as systematics. 

After dividing the lightcurves in seven segments 
of 500 datapoints each, each covering about 10 d 
with a uniform cadence of 30 min, we calculate the 
autocorrelation function and determine the posi¬ 
tion of the first non-zero peak, which is a measure¬ 
ment of the fundamental period in the dataset. 
For 10 objects, we find a consistent period in all 
7 segments, for 5 more the same period is found 
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Fig. 1.— Processing of the SAP lightcurve for EPIC204614722, see text for details. 


for 5 or 6 segments. The scatter in the individual 
period measurements is 1-8%, which we adopt as 
the uncertainty. One additional object shows an 
obvious ~ 5 d period, but also strong additional 
variations, which could be intrinsic to the object, 
but might also be partly instrumental. The best 
periods and the uncertainties are listed in Table 
[2 The lightcurves are plotted in phase to the best 
period (for the first 10 days of the campaign) in 

Fig. 13 

To check for false positive signals due to con¬ 
tamination from other sources, we extracted the 
light curves of K2 stars near our targets as listed 
in MAST (typically much brigher, and within 
~ 2 — 3"). We also tested various apertures on 
a target-by-target basis - from the central pixels, 
from the pixels surrounding them, and from any 
other conspicuously bright pixel in their respective 
postage stamp images - and visually examined the 
extracted light curves. Our analysis indicates that 
the flux modulations in our targets’ light curves 
originate in their respective central pixels only. 


^We note that we find a period of 0.48 d for three tar¬ 
gets not in Table [2 EPIC204367193, EPIC204418005, and 
EPIC204086791. We discard these because the lightcurves, 
phase plots and autocorrelation functions do not look con¬ 
vincing. These could be residuals of the systematics and 
have to be re-examined when new software is available. 


In addition, with the exception of Mars moving 
across the detector containing our targets for a 
few days in October of 2014, there are no obvi¬ 
ous false positive sources. We also confirmed that 
the periods estimated here are also visible in the 
lightcurves made p ublicly available by|Vanderburg| 


& Johnson ( 2014| ^ Therefore we are convinced 


that the periods in Table are of astrophysical 
origin. 

The measured periods range from 0.2 to 5d, 
with a median of 1.1 d. The amplitudes of the 
lightcurves are typically a few percent, with the 
shape usually well approximated by a sinus curve. 
According to spectroscopy analysed in the litera¬ 
ture, the objects with measured periods have spec¬ 
tral types of M5 to M7.5 and effective tempera¬ 


tures from 2500 to 3000 K (Slesnick et al. 2008 


Lodieu et al. 20111, corresponding to masses of 


0.02 to 0.09 M0 (Baraffe et al. 20151. The ma¬ 


jority of these objects are likely at the high end 
of this mass range. Given that our targets are 
all mid to late M dwarfs and thus magnetically 
active objects, the most straightforward explana¬ 
tion for the periodicities are magnetically induced 
cool spots on the surface, which modulate the flux 
over a rotational cycle. Therefore, we interpret the 


^https://www.cfa.harvard.edu/ avanderb/k2.html 
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Fig. 2.— Phased lightcurves for the first 
periods, in the same order as in Table 


ten days of the K2 campaign 2 


for the 16 objects with 


detected 
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Table 1: UpSco rotation periods. 


EPIC 

2MASS 

Kmag 

TeS (K)^ 

P (d) 

Perr (d) 

203420271 

16033799-2611544 

15.167 


1.30 

0.01 

204204414 

16153648-2315175 

16.417 

3025 

0.61 

0.02 

204451272 

16082229-2217029 

15.556 

3025 

0.96 

0.02 

204099713 

16112630-2340059 

15.971 

2818, 3025 

1.80 

0.03 

204344180 

16143287-2242133 

16.775 

2630, 2935 

1.79 

0.06’"’^ 

204614722 

16095217-2136277 

13.933 

2570 

0.66 

0.03 

204105498 

16124692-2338408 

14.964 

2754, 2910 

0.28 

0.02 

204655153 

16105499-2126139 

15.316 

2754 

0.51 

0.01 

204126288 

16164539-2333413 

16.293 

2884, 3025 

0.21 

0.01 

204393705 

16132665-2230348 

16.123 

2910 

1.53 

0.03 

204239143 

16113837-2307072 

15.143 

2910 

0.73 

0.02 

204149252 

16133476-2328156 

14.843 

3025 

0.98 

0.02 

204078097 

16095852-2345186 

13.992 

2935 

1.42 

0.04’^’'= 

204078986 

16080745-2345055 

15.786 

2935 

1.72 

0.08’^ 

204439854 

16113470-2219442 

14.662 

2035 

1.21 

0.05 

203348744 

16030235-2626163 

14.614 


~ 5 

0.5=^ 


“mid-infrared excess emission 

^irregular eclipses in the first part of the lightcurve 
“rapid amplitude changes 


^references for Tefj: Slesnick et al. (2008i, 


Lodieu et al. 


120111 


periods as rotation periods. For a more detailed 
discussion on the causes of periodic variability in 


very low mass objects, see Scholz & Eisldffel (2004 


20051. 


For the other 35 objects in the sample we can¬ 
not determine rotation rates. It is conceivable 
that the periods measured here are biased, partic¬ 
ularly since we are using an activity phenomenon 
(spots) to derive periods, and activity may depend 
on rotation. Ultimately this can only be tested 
by complementary, activity-independent measure¬ 
ments of rotation, such as spectroscopic v sin i. At 
this point there is no evidence for a bias in the 
photometric period samples; the rotation-activity 
relation is flat in mid M dwarfs and so far v sin i 
have confirmed the fast rotation of very low mass 
objects (e.g. Jackson fc Jeffrie^|2012 1. Therefore, 
for the remainder of the paper we will assume that 
these periods are representative of young BDs in 
UpSco. 


The periods are persistent throughout the cam¬ 
paign, i.e. over at least 72 d, with stable ampli¬ 
tudes. Most objects show clear signs of phaseshift: 
If we plot the datapoints as function of phase for 
each of the segments, the maximum shifts through 


the observing campaign by 0.1-0.5 cycles over the 
course of 72 d. Comparing the phased lightcurve 
for various segments also shows in some cases in¬ 
dications for changes in the shape. Thus, while 
the surfaces of our targets are continuously cov¬ 
ered by spots, the parts of the surfaces that induce 
the variability change over time. This is evidence 
for evolving spot pattern on young BDs. In one 
case (EPIC204614722) the period significantly de¬ 
creases over the observing campaign, from 0.69 d 
to 0.61 d, most likely because the variability is in¬ 
duced by spots located at varying and differen¬ 
tially rotating latitudes. 


Out of the 16 objects, 4 have mid-infrared 
colour excess indicating the presence of a disk 
(Class II, Dawson et al.|[20T3 l; the remainder are 
disk-less based on their broadband colours (Class 
III). One of the Class II sources, EPIC204344180, 
shows a series of eclipses in the first part of its 
lightcurve (Fig. left panel). The width and 
depth of the features vary slightly, and the eclipses 
disappear after about 20 d. The dominant period¬ 
icity, however, persists over the entire lightcurve. 
This object might belong to the object class of 
AA-Tau analogs, recently termed ’dippers’, which 
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show either regular or irregular eclipses, presum¬ 
ably caused by a wall at the inner edge of the 


ods for the cluster around e Ori (Scholz & Eisloffel 


disk (Morales-Calderon et al.|[2011 Stauffer et al. 


20151. To our knowledge, this is the first exam¬ 


ple of a BD ’dipper’. The fact that the rotation 
period of this object is consistent with the period 
of the eclipses indicates that the disk feature has 
to be near co-rotation radius. Another one of the 
sources with disk, EPIC204078097, has only in¬ 
significant colour excess at 3-5 ^m, but clearly ex¬ 
ceeds the photospheric emission at 12 ^m, which 
makes this a BD transition disk. This object shows 
irregular amplitude changes throughout the ob¬ 
serving campaign (Fig. right panel), possibly 
due to short accretion bursts, maybe analogous to 
the ones reported for stars by |Stauffer et al. (2014) 
and Findeisen et al. (20131. As pointed out above. 


EPIC203348744, another Class II object, features 
irregular variability which could at least partially 
be caused by accretion or the presence of a disk. 

4. Discussion 

The new periods for UpSco BDs give us an op¬ 
portunity to probe the rotational evolution in the 
substellar regime. At ages of 1-20 Myr, stars spin 
up when angular momentum is conserved, because 
they contract towards their final main-sequence 
radii. Angular momentum losses occur either due 
to magnetic star-disk coupling or accretion pow¬ 
ered stellar winds (Matt et al. 2010, 2012). For 


a review of these processes and the observational 


evidence for disk braking, see Bouvier et al. (2014) 


In Fig. we compare the rotation periods in 
UpSco presented in this paper (Table with the 
two other sizable samples of BD periods, in the 
ONC (Rodriguez-Ledesma et al. |2009| ) and a orjg 
and with simple rotational evolution tracks. The 
models are calculated based on the evolutionary 


models by Baraffe et al. (2015), for a fiducial BD 


mass of 0.07 Mq. Adopting a different mass (0.04- 
0.07 Mq) shifts the tracks slightly. All tracks as¬ 
sume angular momentum conservation. We also 
overplot the physical barrier for the rotation pe¬ 
riod (’breakup period’). 

For tj Ori, we combine the samples by Scholz fc| 
Eisloffel (2004) and Cody & Hillenbrand (2010). 


We also include in the a Ori sample the BD peri- 


■’Joergens et al. 
BDs in Chamae 


I 2003 1 published 3 more periods for young 


2005), a poorly characterised region, but note that 


their inclusion does not change any of the following 
results. For the samples by |Rodrfguez-Ledesma| 


et al. (2009) and Cody & Hillenbrand (2010) we 


adopt a magnitude cutoff of / = 16.5 to separate 
BDs from stars; this is subject to uncertainties 
in evolutionary tracks and cluster distances. For 
the other two samples, masses have been deter¬ 
mined by comparing the IJHK photometry with 
evolutionary tracks in the original papers. In to¬ 
tal we use 179 periods in the ONC and 29 periods 
in tj Ori. We overplot the median and the 10% and 
90% percentile. All samples should be sensitive to 
periods ranging from ~ 0.2 to several days, in the 
case of the K2 sample up to several weeks. 

The choice for the ages of crOri and UpSco is 
relevant for the discussion. For a Ori members the 
plausible age range is 3-6Myr (Sherry et al.||2008 


Bell et al.|2013 ), for UpSco 5-12Myr (Preibisch et 


al.|2002 ' |Pecaut et al. Both regions are sig¬ 


nificantly older than the ONC, and UpSco is prob¬ 
ably again significantly older than tr Ori. In Fig. 
1^ we choose ages of 4 and 8 Myr for cr Ori and Up¬ 
Sco, respectively. Age spread is probably present 
in both regions, but since ages for individual BDs 
are not easily determined, this is neglected here. 

The median and percentiles for the UpSco pe¬ 
riods are slightly lower than in the ONC, consis¬ 
tent with spinup due to contraction. Compared 
with tj Ori, the median and lower limit are actu¬ 
ally slightly higher in UpSco. We do not put too 
much emphasis on this finding at this point, since 
we cannot be sure that we are not missing some of 
the ultrashort periods in the K2 lightcurves, which 
may have been removed together with systematics. 
The upper limit in UpSco is lower than in cr Ori, 
mostly because a Ori has a larger fraction of slow 
rotators with P > 2d (8/29 vs. 1/16). Some of 
the periods in a Ori and e Ori are below or very 
close to the breakup period, this may have conse¬ 
quences on the rotational evolution, see discussion 


Scholz & Eisloffel (2005). 


eon-I (age 1-2 Myr), with P = 2.2 ... 3.4 d. 


From Fig. it is evident that the current pe¬ 
riod census does not require any disk braking to 
explain the evolution of the period median and up¬ 
per/lower limits. Only ages at the upper end of 
the plausible range (6 Myr for crOri, 11 Myr for 
UpSco) may require to include some disk lock¬ 
ing in the evolutionary tracks, but only for the 
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Fig. 3.— Partial lightcurves for two objects with periodic variability to highlight the irregular eclipses 
(EPIC204344180, left panel) and bursts (EPIC204078097, right panel). 



Fig. 4.— Rotation periods for brown dwarfs as a 
function of age, compared with evolutionary tracks 
assuming angular momentum conservation. The 
breakup period is plotted as well. The periods at 
1 Myr (ONC) are randomly spread out in age for 
clarity. 


slow rotators and only over a locking timescale of 
< 2 — 3 Myr. This is about about half as long 
as the locking timescale for slowly rotating low- 
mass stars ( Herbst h Mundt||2005 Gallet & Bou- 


vier 


2013). Thus, independent of the choice of 


the ages, this analysis robustly demonstrates that 
disk braking is inefficient in the BD regime. This 
is consistent with earlier findings of ’moderate an¬ 
gular momentum loss’ in very low mass stars, by 


Lamm et al. (2005). It is notable that the four ob¬ 


jects with disk in the UpSco sample are among the 
six slowest rotators, with periods of 1.4, 1.7, 1.8 d, 
and 5 d well above the median. This indicates that 
disk braking does have some effect, but only for 
very few objects and/or on short timescales. 


5. Summary 

We have measured a sample of rotation peri¬ 
ods for young brown dwarfs in the UpSco region, 
using lightcurves from campaign 2 of the K2 mis¬ 
sion. Our periods range from a few hours up to 5 d, 
with a median just above one day, confirming that 
BDs, except at the youngest ages, are fast rotators. 
Four of the slowest rotators have mid-infrared ex¬ 
cess emission due to the presence of a disk; at least 
two of them show signs of disk eclipses and ac¬ 
cretion in the lightcurves. We compare the new 
periods with previously published samples in the 
ONC and cr Ori and show that the period evo¬ 
lution from 1 to 10 Myr is consistent with no or 
little rotational braking, in contrast to low-mass 
stars. This confirms that disk braking, while still 
at work, is inefficient in the BD regime. This main 
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finding is potentially an important constraint on 
the mass dependence of the braking mechanism. 
Compared with low-mass stars, young BDs have 


much lower accretion rates (e.g. Natta et al. 


weak magnetic fields (Reiners et al. 2009), and 


20041, 


possibly lower ionisation rates at the inner edge 
of the gas disk. These characteristics may influ¬ 
ence the efficiency of rotational braking. 


We thank Thomas Barclay for his assistance 
with PyKE and Andrew Vanderburg for a fruit¬ 
ful discussion on the details of his data reduction 
technique. This work was supported in part by 
NSERC grants to R. Jayawardhana. This work 
made use of PyKE ( Still fc Barclay|[2012 ), a soft¬ 
ware package for the reduction and analysis of Ke¬ 
pler data. This open source software project is 
developed and distributed by the NASA Kepler 
Guest Observer Office. 
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